Quantum correlations in inflationary spectra and violation of Bell inequalities 



David Campo^ and Renaud Parentani^ 

^Department of Applied Mathematics, University of Waterloo, 
University Avenue, Waterloo, Ontario, N2L 3G1 Canada 
^Laboratoire de Physique Theonque, CNRS UMR 8627, 
Bdtiment 210, Universite Paris XI, 91405 Orsay Cedex, France 

In spite of the macroscopic character of the fluctuation amplitudes, we show that the standard 
inflationary distribution of primordial density fluctuations still exhibits inherently quantum me- 
chanical correlations (which cannot be mimicked by any classical stochastic ensemble). To this 
end, we propose a Gedanken experiment for which certain Bell inequalities are violated. We also 
compute the effect of decoherence and show that the violation persists provided that the decoher- 
ence lies below a certain non-vanishing threshold. Moreover, there exists a higher threshold above 
which no violation of any Bell inequalities can occur, so that the corresponding distributions can be 
interpreted as stochastic ensembles of classical fluctuations. 



Inflation implies that the primordial density fluctu- 
ations arise from the amplification of vacuum fluctua- 
tions As a result, the initial vacuum state becomes 
a product of highly squeezed two-mode states char- 
acterized by the wave vectors k and — k. This follows 
from the homogeneity of the state and the linearity of 
the mode equation. However the linear treatment is 
only approximate. Hence, when no longer neglecting 
the weak non-linearities of inflationary models 3] , and 
when restricting attention to a given two-mode sec- 
tor, decoherence will unavoidably develop. Since the 
interactions preserve the homogeneity and are weak 
(so as to justify the use of a Hartree approximation) 
inflationary distributions belong to the class of Gaus- 
sian homogeneous distributions obtained by slightly 
decohering the standard distribution obtained by ne- 
glecting non-linearities 4, 5, 6]. 

Given the macroscopic character of the fluctua- 
tion amplitudes at the end of inflation, one might 
ask whether some quantum properties have been pre- 
served by these distributions, and if so, how to distin- 
guish those which have kept some from those which 
have lost them. To show that the standard distri- 
bution has kept quantum properties and to sort out 
these two classes of distributions, we shall analyze the 
possibility of violating Bell inequalities 0- 



I. THE INFLATIONARY SPECTRA 

A. The 'standard' distribution 

In inflationary models based on one inflaton field, 
one can express the dynamics of the linear pertur- 
bations around the background solution of Einstein 
equations in terms of a massless scalar field minimally 
coupled to gravity, 0. This field is the gauge-invariant 
linear combination of the gravitational potential and 
the infiaton fluctuations which satisfies canonical com- 
mutation relations We briefly outline how one 
obtains the two-mode squeezed states and the time- 



coherence of the modes 0, 1^1 ■ 

The background geometry is a Friedman- 
Robertson- Walker space-time, with line element 



(1) 



The spatial sections in the present Hubble volume are 
flat to an excellent precision if inflation has lasted 
enough to have increased the radius of curvature to 
super-horizon scales. The scale factor, as a function 
of the conformal time rj is given by a{ri) = —1/Hrj 
in the inflationary period we approximate by a de 
Sitter space. At the reheating rjr, we paste this 
law to a radiation dominated phase characterized by 
a{rf) = {f] — 2r]r)/ {Hrjf). We can then follow the time 
dependence of modes from the onset of inflation till 
horizon re-entry in the radiation dominated era. 

After expanding the action of matter and metric 
fluctuations up to quadratic order around the back- 
ground solution, the Hamiltonian for each complex 
Fourier mode (/)k is 



\drj(t)U 




(2) 



The state of (py^ at the onset of infiation follows from 
the kinematics of the background ;9j: if some inho- 
mogeneities are present at an early stage of inflation, 
they are diluted by the quasi-exponential expansion 
over proper lengths bigger than our causal patch, i.e., 
to us, they are indistinguishable from the background. 
One is thus left with the vacuum. One can also reach 
this conclusion by considering the evolution of pertur- 
bations backward in time. Then, non-vacuum contri- 
butions (classical inhomogeneities as well as quantum 
excitations of the fleld) to the energy-momentum ten- 
sor blow-up thereby violating the consistency condi- 
tion of no-large backreaction. In conclusion, all the 
modes observable today were in their ground state 
about 70 e-folds before the end of inflation (the min- 
imal duration of inflation to include today's horizon 
inside a causal patch). 
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This ground state is unambiguously defined for the 
relevant modes because, at that time, their wave- 
length was much shorter than the Hubble radius (i.e. 

3> df^a/a in Eq. ((2Jl). One therefore works with 
the vacuum defined in the asymptotic past by posi- 
tive frequency modes 

lim [za,-A;]C(77)-0, (3) 

77-^ — 00 

where (f)'^^{r]) are solution of the mode equation: 

9^0k + (^k^ - <^k - , 

Mv) = aJT^r iv) + a!_"j0r * iv) ■ (4) 

The time-dependent term d^^a/a mixes the positive 
k and negative frequency — k modes. Therefore the 
vacuum unitarily evolves into an entangled state: a 
two-mode squeezed state. Confusion about the origin 
of this entanglement seems to persist: entanglement 
is not a consequence of the reality (self-adjointness) 
of 4i{ri,:ic). Pair production of charged particles (in 
an electric field) also leads to entangled states. Two- 
mode entanglement is due to the mixing of positive 
and negative frequencies and to the homogeneity of 
the background. 

At the end of inflation, (?^a vanishes and the ampli- 
fication stops. During the radiation dominated era, in 
the absence of iso-curvature modes, the fluctuations 
of physical fields (those of the gravitational potential 
and of cold dark matter and radiation density fluc- 
tuations) are all linearity related to the values of </> 
and dr/cl) evaluated at the end of inflation. This is true 
both classically and quantum mechanically because of 
the linearity of the transfer matrix. 

Instead of working with these physical fields, it is 
simpler to keep working with (jik, since its time evo- 
lution is simply governed by e^*'^''. We shall use this 
behavior to further simplify the writings by introduc- 
ing creation and destruction operators fik, aj^ defined 
in the adiabatic era by 0k = (cke"*'^'' -|- h.c.)/V2k. 
In terms of these operators, the 'in' vacuum (i.e. the 
state annihilated by the aj^" of Eq. Q) reads 

k 

X exp (zfc 4 alk) |0' ^) ^ |0' "k) • (5) 

The tilde tensorial product takes into account only 
half the modes. It must be introduced because the 
(squeezing) operator exp (^Zk a^. aLk) ^''^^s on pairsl 

of modes. The complex number Zk characterizes the 
two-mode states of wavenumber k — |k|. It is related 
to the following expectation values 

(4ak') =nfc5'(k-k') , (akflk') =Cfc,53(k + k').(6) 



The mean occupation number is = /(|zfe| — 1) 
and the cross correlation term is Ck — Zk/{\zk\ — !)• 
They obey \ck\ = nk{nk + 1). This is the highest 
value of |cfc I for a given n^. It is obtained for the states 
which minimize Heisenberg uncertainty relations @. 

Asking that the spectrum describes fluctuations 
with r.m.s. amplitude 10^^ for wavelength observable 
in the CMB, one gets the rough estimation Uk ^ 10^"°. 
The phase of Zk determines the temporal phase of 
the modes as they re-enter the Hubble radius |^. 
Namely, the power spectrum of 0k at conformal time 
r] is given by a product of two modes proportional to 
gifc?) _|_2^g-»fc'7_ jn inflationary models, Zk is very close 
to —1, giving rise to sin(fc77) functions, i.e. to 'growing 
modes'. 



B. Partially decohered Gaussian distributions 

The two-mode entanglement of squeezed states fol- 
lows from the homogeneity of the vacuum and the 
linearity of the mode equation. However, as discussed 
in the introduction, the linear treatment is only ap- 
proximate. When no longer neglecting the weak non- 
linearities of gravitation |3| , and when restricting at- 
tention to a given two-mode sector, some decoherence 
develops. This is very similar to the description of 
an experiment of diffraction of a slightly anharmonic 
cristal. The calculation of the reduced density ma- 
trix of the two modes under scrutiny is very diffi- 
cult. But the weakness of the interactions guaran- 
tees that the modifications of the power spectrum are 
extremely small. It also guarantees that the reduced 
density matrix can still be accurately described by a 
Gaussian distribution. Moreover, since homogeneity 
is preserved, these distributions always factorize into 
products of uncorrelated two- mode sectors. In conclu- 
sion, inflationary distributions of primordial fluctua- 
tions belong to the class of Gaussian homogeneous dis- 
tributions obtained by (slightly) decohering the stan- 
dard distribution pin = |0, m)(0, in\ 

These are described by products of two-mode den- 
sity matrices p2 which are characterized by the ex- 
pectation values where nj, and \ck\ are now in- 
dependent parameters. We shall work at flxed nk 
(fixed power spectrum), and fixed arg(cfe) (changing 
this phase does not change the entropy p]), and let 
the norm \ck\ vary. A convenient parametrization is 
provided by Sk- 

|cfc|^ = (rifc + l)(nfe - Jfc) . (7) 

Squeezed states ^ are maximally coherent and cor- 
respond to Sk — 0. The least coherent distribution, a 
product of two random density matrices, corresponds 
to Sk = nk- Hence, the parameter Sk controls the level 
of decoherence of the distribution, or equivalently the 
strength of the correlations between the k and the — k 
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sectors. It also fixes the r.m.s. amplitude of the de- 
caying mode, as seen by decomposing the field modes 
in terms of the growing and the decaying modes: 



sui{kri) - cos{kr]) 



Vk 



(8) 



When (5 < n, 
(Skffk) 

{dkdl) 



n + i - Re{c) ^2n(^l + 0(^)^ , 
n + i + i?e(c) = n^/-* (^l + 0(^)^ , 
/^(c) = ^- + 0{n^l^) « (gk.gl> . (9) 



We have dropped the indexes k since each two-mode 
can be analyzed separately (unless one is considering 
localized wave packets). 

We emphasize that the increase of entropy from the 
initial pure vacuum state to a statistical mixture fol- 
lows, as usual, from the neglect of correlations, see 
footnote p. 188 in 0. In other words, ignoring the 
correlations between several two-modes is the coarse 
graining we adopt. We also notice that this Gaus- 
sian ansatz is a first order (Hartree) approximation. 
With the future experiments, deviations from Gaus- 
sianity will hopefully be measurable and will become a 
powerful tool to raise the degeneracy between various 
models for the primordial universe 
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II. CLAUSER-HORNE INEQUALITIES FOR 
INFLATIONARY SPECTRA 



A. Local Hidden Variable completions of 
Quantum Mechanics and the Clauser-Horne 
inequality 

Bell inequalities and their generalizations form a 
set of constrains on the statistics of outcomes of mea- 
surements when one demands that these outcomes are 
compatible with the classical concept of locality T?] . 

Consider a system made of two subsystems I and 
II and let the two parts interact. In general they will 
become correlated. Suppose that after 'turning off' 
the interactions, the system is in a pure state of the 
form 



\X) =Y.^^^)n®\^)n ■ 



(10) 



where there is more than one term in the sum, where 
and IV')^^ s-re two sets of orthonormal states in 
the Hilbert spaces and of I and II respectively, 
and where P„ = |A„p < 1 are the probabilities for 
the n-st outcome to be realized such that 'Ylm = 1 • 
The states (IIUII are a special case of entangled states, 
see Section III.C for the definition of entanglement 



for statistical mixtures. (Notice that each two-mode 
squeezed state in Eq. Q is an entangled state.) Being 
orthonormal, the \<^\ and |V')i^ can be seen as eigen- 
states of observables ® 1^^ and \^ ® acting 
non-trivially on Ti,^ and Ti^^ respectively. 

After having let the two subsystems interact, we 
separate them to two distant places where two mea- 
surement separated by a space-like interval will be 
performed. We suppose that we have a ensemble of 
copies of the same system all prepared in state \x) in 
order to perform different measurements and to ac- 
cumulate statistical data. The space-like character of 
the interval guarantees that the state of one subsystem 
cannot be affected by any measurement performed on 
the other subsystem. Once the apparatus I has mea- 
sured the observable on system I, we can predict 
with probability 1 (as checked by reproducing the ex- 
periment in the same conditions a large number of 
times and bringing the data together) the outcome 
of the measurement of 1^ (E) 5*^^, irrespectively of the 
values of the probabilities Pn and despite the space- 
like separation of the events of the measurements on 
subsystems I and II (supposing we have the a priori 
knowledge that the system is in the entangled state 
Ix), as we do in inflationary cosmology, see eq. (6)). 

In order to explain these correlations one might be 
tempted to consider that they result from unknown 
common properties shared by the two subsystems, 
and whose values were assigned while I and II were 
still interacting. These properties, to which we do 
not have access by the experiments, are called Hidden 
Variables. The Hidden Variable Programm consists 
in establishing the conditions which have to be met. 
In particular. Bell's theorem and its generalizations 
0, 0, Q| show that Quantum Mechanics cannot be 
embedded in a Local Hidden Variable Theory. In this 
paper, we will consider a particular type of inequalities 
called the Clauser-Horne (CH) inequalities 
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1. Local Hidden Variable Theory 

A Local Hidden Variable Theory (LHVT) is the 
specification of 

i) the conditional probability p$ [0; A] (resp. A]) 
that the measurement of (resp. ^'^^) gives the 
eigenvalue 4> (resp. tp), given a value A of the hidden 
variables, 

ii) a probability measure p{X)dX. 

The requirement of locality forces that 

1) the conditional probabilities p$ (resp. p^s,) is inde- 
pendent of 4'^^ (resp. $^), 

2) the probability distribution p{X) is independent of 
both and 'I'^^, 

3) given a value A of the hidden variables, the outcome 
of the measurements of $^ and ^'^^ are statistically 
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independent, i.e. 

V'; A) A] X A] . (11) 

A given quantum state p e Ti,^ <E) H^^ admits a 
LHVT if, and only if, for any pair of observables 
with spectral decompositions — J2n 'I'nf^t ^^'^ 
= X^nV^nll*, the statistical predictions of the 
model reproduce those of Quantum Mechanics, i.e. 
the joint probability is 



P<5*(0, V') 



dXp{X)p^^{(l),'ip; A) 
Tr(pfl*flf 



(12) 



2. Clauser-Horne inequalities 



The Clauser-Horne inequality is a condition veri- 
fied by the joint probabilities Pii,^(4>,ij}) if Quantum 
Mechanics can be embeded in a LHVT. They follow 
from the algebraic inequality satisfied by any set of 
four real numbers (x, x' , y, y'), all lying in the interval 
[0, 1], see [Lsl for a proof: 



xy -f xy' -\- x'y — x'y' — x — y < . 



(13) 



Identifying x = p*[V';A],a:' = p^s,['^p■, X], y = 
pqi[ip] X], y' — pqi[ip'; A], and averaging over A one gets 

[P^ * (0, ^A) + P<t * (0>') + -P* * (0' , ^) 
- P$*(0',^')] X [P*(<^)+P*(V')]"' < 1, (14) 

where P$, P-j, are the marginal probabilities. 



We want make a pose and advertise the perti- 
nent role of coherent states in inflationary cosmology. 
Firstly, for quadratic Hamiltonians, they are minimal 
uncertainty Gaussian states which are stable in time. 
They are thus the quantum counterparts of classical 
points in phase space. More generally, one expects 
that they are still stable as long as the secular ef- 
fects of non-linearities (changes in the width of coher- 
ent states, formation of small scale structures -whirls 
and tendrils- of the classical trajectories remain 
small. Such is the case in cosmology as long as the 
the linear regime (and therefore the Gaussian approx- 
imation) are good approximations. 

Secondly, for quadratic systems weakly interacting 
with a reservoir of modes (here the coupling is through 
the non- linear gravitational terms), there exists a fi- 
nite time after which the system reaches a state which 
is statistical mixture of coherent states (more gener- 
ally, of minimal uncertainty Gaussian states) 0, [T9l | 
24]. 

Coherent states can therefore be thought as provid- 
ing a particular realization of the Gaussian ensemble 
describing the primordial fluctuations. The probabil- 
ity that the semi-classical values v and w of the k-th 
mode be realized is 



Q{v, w; S) 



TT[p2{S)n{v,w)] 



1 



71+1 



exp 



(n + 1) 



l + <5 



exp 



\w — w{v)\^ 
l + d 



(16) 



B. A set of observables for Gaussian states 

In searching for a violation of Bell inequalities in in- 
flationary spectra, one encounters an unexpected dif- 
ficulty because we are dealing with Gaussian states. 
Indeed, observables which are polynomes of the field 
amplitude and its conjugate momentum cannot vio- 
late Bell inequalities. The reason is that the Wigner 
representation of Gaussian density matrices is posi- 
tive. Hence, it can be used as a probability distri- 
bution function to build an LHV description of the 
expectation values 0. 

We must therefore consider observables which do 
not have a direct classical counterpart. Nonetheless 
the operators we shall consider have a clear and sim- 
ple meaning. They are projectors on a given pair of 
coherent states |w,k) and jw, — k) 

n(w, w) = |u, k) (v, k| (g) \w, -k) (w, -k| . (15) 

The coherent states obey, by definition, ak|w,k) = 
v\v, k) and a_k|w, — k) = wlw, — k). 



We have written this probability in an asymmetric 
form to make explicit the power spectrum {= n -\- 1), 
and the much smaller width (= I + S) governing the 
dispersion of the values of w around w{v) = v*c/{n -\- 
1). Hence, as long as 5 <C n, the amplitude of the 
mode — k conditional to the detection of the mode k 
in the coherent state \v) is fixed by that v. We see 
explicitly that it is 6, the decoherence level, which 
governs the strength of the correlations between the 
k and the — k sectors. 

Thirdly, given a realization, one can calculate con- 
ditional correlations. The latter have specific spatial 
properties best revealed by wave-packets. Although 
their Fourier content and their localization depends on 
the chosen wave-packet, their space-time structure is 
uniquely determined by the frequency mixing. More 
precisely, given the ensemble of realizations selected 
by Htu.p. the projector on some wave packet, the con- 
ditional field amplitude is (see for more details) 

(0(77, X)) = ^ =<Pr + (Pl, (17) 
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where 



(Pk 
(27r)3/2 

'(Pk 
(27r)3/2 



Ke'''^0r(^) + c.c.) , 
KK)e-*''</.rW + c.c.) 



One sees that the centers of these wave-packets prop- 
agate in opposite directions on the future Ught cone 
from the event of creation at the reheating. Hence, 
on a given time shce, e.g. the Last Scattering Surface 
(in the matter dominated era), they are separated by 
more that twice the Hubble radius at that time. 

The important conclusion for Bell inequalities is 
that the measurement (|15|l (more precisely, its gen- 
eralization which contains a product over k so as to 
describe a localized wave packet) can be performed 
on subsystems which are causally disconnected. Vio- 
lations of Bell inequalities with such observables were 
first considered by Wodkiewicz et al in the context of 
Quantum Optics [l^. 



C. Separable states 

A two-mode state is called separable, or classically 
correlated, if and only if it can be written as a convex 
sum of products of one-mode density matrices |2Q| |. 
Otherwise, the state is called entangled. Restricting 
to the class of homogeneous Gaussian states (O , sep- 
arable states are of the form 

Pr-(S)^ f^^P{v,w;6)Il{v,w), (18) 



where P is a Gaussian function. For homogeneous 
distributions, it is given by (Qj 



P{v,w;6) = 



where IS! = r? — 



X exp 



\w — [cv* ln)Y 
A'/n 



must be positive. This implies 
|c| < n, or (5 > n/(n -I- 1) > 1. In other words, all the 
states such that n <\c\< \/n[n -I- 1) are entangled. 

The physical meaning of separable states is revealed 
by the fact that they can be prepared by a classical 
protocol: a random generator produces the four real 
numbers encoded in (u, w) with probability P. The re- 
sult of each drawn is send by classical communication 
channels to two distant observers performing separate 
measurements on the subsystems k and — k respec- 
tively so as to prepare them into the two-mode coher- 
ent state |u,k)|w,— k). By construction, the statisti- 
cal properties of separable states can be interpreted 
classically. In particular, separable states cannot vio- 
late Bell inequalities as demonstrated in j2Q| . We will 
use them to derive the CH inequalities. 



When applied to separable Gaussian states H18I) , the 
probability H16|l is 

Q-P- {v,w;6)= f—^ P{a, P; 5) p{v\a) p{w\p) , 

J TT TT 



where p{v\a) — |(w|q;)|^ — e^l^""^! < 1. This has 
the structure of Eqs. l(TT|l and l(T^ . where {a, f3) play 
the role of the "hidden variable" A, and P{a, f3; S) the 
role of the probability density p(A). Then Eq. ifTljl 



gives 



Civ,w;S) = [Q{0,0;S) + Q{v,0;S) 



(19) 



III. 



n + 1 

+ Q{0,w;5)-Q{v,w;6)]x ^— <1. 



VIOLATION AS A FUNCTION OF THE 
LEVEL OF COHERENCE 



We now have all the elements to search for combi- 
nations of V and w which maximize C of Eq. H19|l . The 
maximum is reached for aTg(c*vw) = tt and = \w\. 
To make contact with 0|, we choose the arbitrary 
phase 2arg('i;) = arg(c) to get w — ~v. The maxi- 
mum of this function is reached for 



\vMiS)\' 
1 + S 



In 



n—5 
n+1 



1 - 

In 2 
~3~ 



1-5 
i+1 



and its value is 



CAliS) 



1 



2(1 + (5) 



1 + 



3 

2^73 



O 



1 + S 



(20) 



(21) 



There is a violation even for macroscopic occupation 
numbers. This fact demonstrates that n ^ 1 is a 
necessary but not sufficient condition to have a clas- 
sical (stochastic) distribution. The latter are only the 
separable ones for which S > I. 

More precisely, the inequality ((T^ is violated for 



<5< 



-1 + 3/2^/3 



0.095. 



(22) 



irrespectively of the value of n when n ^ 1. Even 
though each probability in (|19|l decreases like 1/n, the 
range of S is asymptotically constant in the regime 
relevant for inflationary cosmology. One also verifies 
that, as expected, the maximal violation is obtained 
for (5 = 0, i.e. for pure states. In Fig. 1 we have plot 
C{v) for three values of S. 
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FIG. 1: The loss of violation as decoherence increases. 
The function is normalized to its maximal classical value: 
C = 1. The variable x = gives the square ampli- 
tude of the field configurations. The occupation number 
is n = too, and the three values of 5 are (red, upper), the 
pure state, 0.05 (blue, middle), and 0.1 (brown, lower), the 
border regime. The function C{v,—v,5) is asymptotically 
independent of n in the large n limit. 



IV. FINAL REMARKS 



First, there might be operators which violate Bell 
inequalities for 0.1 < 6 < 1. However it is not guar- 
anteed that these will have a simple interpretation in 
cosmological terms. 

Second, an observational verification of the viola- 
tion would require several steps. First, one should be 
able to distinguish k from — k configurations, i.e. to 
have access to the primordial velocity field, see 0. 
Second, one should isolate rare realizations, far from 
the r.m.s. values, and specified with a precision given 
by the spread of the coherent states {— 1) which is 
also much smaller than n. Given that n ~ 10^*^°, an 
observational verification seems therefore excluded. 

It is nevertheless conceptually important to real- 
ize that if non-linearities, and hence decoherence, are 
weak enough so that ^ < 1, the distribution of primor- 
dial density fluctuations would have kept its quantum 
properties in spite of the macroscopic character of the 
amplitudes and even though we cannot observation- 
ally verify it. 
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